High energy scales in electronic self-energy imaged by optics 
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We use a new technique to directly extract an estimate of the quasiparticle self-energy from the 
optical conductivity which can be easily related to both theory and angle-resolved photoemission 
spectroscopy (ARPES) experiments. In the high T c cuprate Bi-2212 we find evidence for a new high 
energy scale at 900 meV in addition to the two previously well known ones at roughly 50 and 400 
meV. The intermediate scale at 400 meV has been recently seen in ARPES as a large kink which 
optics finds to be weaker and shifted. In YBCO, the three energy scales are shifted to lower energy 
relative to Bi-2212 and we observe the emergence of a possible fourth high energy feature at 600 
meV. 
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One of the most beautiful aspects of studying con- 
densed matter systems is that a large number of very 
different types of experiments are available to probe the 
many-body state, and in the quest for microscopic un- 
derstanding all of the various probes must concur about 
the features observed. This provides a stringent test of 
theories and further allows for the perfection of experi- 
mental technique and analysis. Consequently, if there is 
an opportunity to bring two different techniques closer 
together for the purpose of comparison, both theory and 
experiment will benefit and there will be a great enhance- 
ment of the ability to identify potentially fundamental 
and robust signatures of the underlying physics. Two 
powerful techniques which have been brought to bear 
on the problem of high T c superconductivity are angle- 
resolved photoemission spectroscopy (ARPES) 0,0 and 
far-infrared optical spectroscopy These techniques 
have blossomed with the demands of this area of study re- 
sulting in unprecedented resolution and new methods for 
extraction and analysis of information, such as through 
inversion techniques. 0, IE IE IE IE E| 

Applying inversion techniques to the optical conduc- 
tivity has resulted in the identification of a coupling of 
the charge carriers to a resonance mode at about 50 
meV and a spin fluctuation background extending be- 
yond 400 meV.IEIEIIEGil ARPES, on the other hand, 
has mainly confirmed the low energy scale as, until re- 
cently, they were limited to presenting data for the quasi- 
particle dispersion curve up to about 200-300 meV in 
energy. [IlHIIllIEIIElHlIEIIllIl However, new 
advances have extended this r ang e to greater than 1 eV 
and several works 0, |^E HE 0, HE] are reporting fea- 
tures at high energy at about 400 meV and 800 meV, 
which are speculated to arise from exotic physics. Con- 
firming the existence of these new energy scales and in- 
vestigating a wide range of materials is essential to the 
development of these or other new ideas. 

In this letter, we present complementary results ob- 



tained from optical data and demonstrate, at the same 
time, a new technique which allows for the direct, al- 
though approximate, extraction of the quasiparticle self- 
energy from the optical conductivity, which can be com- 
pared most easily with the ARPES via a dispersion-type 
curve. This is important because optics is a bulk probe 
as opposed to ARPES which is only sensitive to the sur- 
face. In addition, optical spectroscopy measures the self- 
energy due to correlations directly while ARPES finds, 
instead, the renormalized energy and the self-energy fol- 
lows only once the bare dispersion is known. Using our 
new technique to analyze the data on Bi-2212 samples of 
varying doping, we can show that optics also sees features 
at high energy and supports some of the observations seen 
in ARPES, but not all. Moreover, as optical experiments 
are more flexible in the type of materials which can be 
studied, our work goes beyond Bi-2212 to examine the 
high energy scales in YBCO. 

We begin by discussing the theoretical basis for the 
new analysis of optics which allows for the extraction 
of the quasiparticle self-energy Within the experimen- 
tal community, it has become common to analyze the 
complex optical conductivity cr(u>) for the correlated 
charge carriers in terms of the optical self-energy, E op = 
Y±° p + iT,2 P {oj) , using an extended Drude model of the 
form: 0, [13 
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where uj p is the plasma frequency. S op can be written 
in terms of the optical scattering rate l/r op (w) and the 
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optical mass renormalization 1 + \ op (u>) by 
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and it is routine to present quantities such as 1/t° p or 
S op . Despite the suggestive structure of this form, E op is 
not the quasiparticle self-energy S 9P and 1 /t° p is not the 
quasiparticle scattering rate. However, these quantities 
do encode such information and within a set of reasonable 
approximations, a surprisingly simple procedure can be 
applied to extract Tfl p from E op . 

If we can assume, for the sake of argument, that the 
renormalization and scattering of charge carriers is due 
to some sort of frequency-dependent interaction with a 
bosonic spectrum which we will call a 2 F(uj) in analogy to 
the conventional electron-phonon spectral function, then 
we can show that El IHElEil 
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where the electronic density of states N(oj) can be 
energy-dependent and we have ignored the difference be- 
tween the transport and quasiparticle spectral densities, 
a 2 r F and a 2 F, respectively. The approximations behind 
these formulas and their range of validity may be found 
in the work of Allen [30|. Shulga et al.|3l|. and Sharapov 
and Carbottejs^- For convenience, these equations are 
stated here for T = and for no impurity scattering. In 
the cuprates, we can assume the clean limit. Finite tem- 
peratures effects arc easily incorporated when needed. If 
we define a new quantity, following the work of Carbotte 
et al.|^, which we will call the optically-derived quasi- 
particle self-energy, 



Y, op - qp (u) = —\u,E op (lo)} 
duj 



(7) 



then from Eqs. JSJ) and iJBJ, the derivative in Eq. Q gives 
exactly Y, qp . 

To illustrate the effectiveness of Eq. {7J to determine 
the exact quasiparticle self-energy, we show a model cal- 
culation in Fig. ^ where an a 2 F(uj) spectrum consisting 
of three delta functions has been used to calculate both 
the exact T, qp from many-body theory and the optical 
conductivity in the usual Kubo approximation. [2ft l27l | 
The optical conductivity is then processed via Eqs. 
® and to obtain Y, op ~ qp . In the figure, we have 
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FIG. 1: (color online) Theoretical illustration based on the 
model described in the text which demonstrates that the 
quasiparticle dispersion curve can be extracted from optics 
(op-qp) via our new procedure. It agrees well in magnitude 
and detailed structure with the exact theoretical curve (qp). 
If the optical self-energy is used, the resulting dispersion curve 
(op) is structureless and does not agree with the exact result. 
The inset illustrates for a real material Bi-2212 the difference 
between the dispersion curve derived from the optical self- 
energy and that from the new procedure which gives rise to 
detailed structure and specific energy scales. 



chosen to plot the real part of £ in terms of a renormal- 
ized dispersion curve where Ek — £fe — Si(-Bk) [34| as we 
wish eventually to make comparison with ARPES in the 
same format. The dispersion curve for Yfl p shows very 
clearly three structures due to the three J-function peaks 
in the a 2 F at 0.04, 0.09, and 0.19 eV 27}. However, the 
curve formed from E op is fairly featureless and does not 
capture any of the underlying structure of the spectral 
function. In addition, this calculation has included a fi- 
nite band density of states jl^, 113, which causes the 
renormalized dispersion curve to cross the bare disper- 
sion curve (a fact that is not well-known, but is a signa- 
ture of finite bands) and once again the curve from S op 
has not crossed the bare dispersion at the same energy 
and indeed the crossing is higher in energy by a factor 
of roughly two. On the other hand, the optically-derived 
quasiparticle self-energy Y, op ~i p gives a dispersion curve 
in excellent agreement with the exact one and it repro- 
duces all of the structures at the correct energies and 
with the proper magnitude. Aside from a small over- 
shoot of the curve between about 0.3 and 0.4 eV due to 
the derivative technique which is known to produce such 
characteristics just after the end of the boson spectrum, 
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FIG. 2: (color online) (a) Dispersion curve derived from optics 
for optimally-doped Bi-2212, for T = 100 K, in a color map 
showing the relative height and width of A(k, w), as is done in 
ARPES. (b) Dispersion curve from (a) repeated (op-qp) and 
compared with the result reported by ARPES and with the 
theoretical LDA curve from Ref. [36| . The grey shading gives 
the width of the optically-derived quasiparticle curve based 
on the imaginary part of the self-energy. 

the op-qp curve also crosses the bare dispersion at the 
same point as the qp curve. The agreement is remarkable 
and similar agreement is found in the detailed compar- 
ison of the real and imaginary parts of Tfl p and Yj op ~ qp 
(not shown). The grey shading shows the broadening of 
the dispersion curve which occurs with E2 and is indi- 
cated by ± I Sip I- The inset of this figure shows the pro- 
cedure of Eq. @ applied to experimental optical data for 
optimally-doped Bi-2212 at T = WOK [TflU^. Again, 
Y°p gives a featureless dispersion curve, but T, op ~ qp re- 
stores structures in the curve (such as the kink at ~ 50 
meV) which can now be compared to ARPES. Fig. 
emphasizes that it is "E op ~ qp which must be used for this 
purpose. A point that must be made before continu- 
ing to such a comparison is that optics is a transport 
property and it does not contain momentum-dependent 
information. Consequently, along with the assumption 
that a 2 r F — a 2 F, we stress that the "dispersion curves" 
presented here are to be viewed as momentum-averaged 
curves in terms of the self-energy and the only remaining 
momentum-dependence enters through the bare disper- 
sion efe which is taken from LDA calculations. 

Having established the technique, we now present our 
new results on high energy scales in Bi-2212 and YBCO 
as imaged by optics. In Fig. [21 we present the analy- 
sis of the optical data for optimally-doped Bi-2212 at 
100 K. The derivative needed in Eq. J7J) to extract 
Ypv-qv can introduce noise and we have applied smooth- 
ing techniques as is reasonable since we are interested 
mainly in capturing global features. For comparison 
with ARPES, we have presented the T, op ~ qp in the form 
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FIG. 3: (color online) Upper frame: The real part of the 
optically-derived quasiparticle self-energy for three Bi-2212 
samples, overdoped, optimally-doped and underdoped, and 
one underdoped YBCO in the ortholl phase [38|. all at T = 
100 K. Lower frame: The imaginary part of the optically- 
derived quasiparticle self-energy. Various energy scales are 
indicated on the curves and discussed in the text. The inset 
gives the a 2 F(u>) obtained in the optimally doped Bi-2212 
from inversion of optical data. [Til ] 

of a dispersion curve using an LDA calculation for 
the bare dispersion e^. With this information, we can 
form a color map of the quasiparticle spectral function 
A(k,u) = - e k - Yf l p ~ qp ) 2 + (£°^) 2 ] 

to connect to the raw data as presented by ARPES. In 
Fig. [2b, we compare our data with that of Graf et al.|2lj 
from ARPES in the nodal direction on a similar optimally 
doped sample. In the ARPES, they identify three energy 
scales with these data at Eq ~ 50 meV, E\ ~ 0.4 eV, and 
E2 ~ 0.8 eV. In our data, we can confirm a structure at 
E2 ~ 0.9 eV and a much stronger kink at Eq ~ 50 meV, 
however, we do not find evidence for a very strong fea- 
ture at ~ 0.4 eV, although the Y, op ~ qp does show a local 
maximum around 0.25 eV. The broadening that we find 
is within the range that is found by ARPES. Thus, we 
conclude that overall, we find agreement with ARPES 
using our analysis as shown in Fig. [2b, however, there 
remain some discrepancies which we will now address. 

Inversion^, |jj of optical data to recover the electron- 
boson spectral density of Eqs. JSJ and © shows that 
the kink at ~ 50 meV may be attributed to a peak in 
a 2 F(ui) which is followed by a broad spin fluctuation 
background extending beyond 0.4 eV (see inset of Fig. [3}. 
This spectrum accounts for the deviations, extending up 
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to ~ 0.6 eV, of the op-qp data from the LDA curve seen 
in Fig. |5Jd. The renormalization found in the ARPES 
data, while larger than that found in the op-qp curve 
around 0.4 eV, is smaller below the 50 meV kink. The 
lack of quantitative agreement between the optics and 
ARPES in these two regions may reflect (i) the differ- 
ences between the quasiparticle and transport spectral 
densities, not included in this work, (ii) the differences 
in k-dependence between optical data (averaged over the 
Fermi surface) and that of ARPES (momentum specific), 
and (iii) that formally the self-energy S(k, ui) can depend 
on the magnitude of k as well as direction, but this goes 
beyond the theory used here. It is clear that renormaliza- 
tions are present in both experiments on a broad energy 
scale up to the end of the spin fluctuation spectrum (~ 
0.4 eV) which we take here to be the E\ scale rather 
than the position of the peak in the self-energy. Further- 
more, we associate the higher energy scale around 0.8- 
0.9 eV with the commencement of a new higher energy 
band (see discussion below). One of the results, strongly 
highlighted in some of the ARPES work[H I2I lisj] is 
the vertical dispersion seen in the data from E\ to Ei- 
The question arises as to whether this is due to some 
sort of ordering^], |3| or if it is a renormalization ef- 
fect. Indeed a vertical dispersion can be obtained within 
standard electron-boson models jH,E3 with appropriate 
choice of parameters and in this instance, the high energy 
scale of this feature points to a spin-fluctuation spec- 
trum rather than a low energy boson spectrum. Thus, 
we would conclude that the large kink is due to renormal- 
ization rather than something more exotic, as was also 
concluded by Valla et al.[23^. 

In the upper frame of Fig- El we show Yf p ~ qv for three 
Bi-22I2 materials: underdoped (UD69), optimally-doped 
(OPT96) and overdoped (OD60), and one underdoped 
YBCO material (UD59) 0, all at T = 100 K in the 
normal state. For the Bi-2212 samples we have indi- 
cated possible features that would correlate with the Eq, 
Ei , and E2 seen in ARPES providing possible confirma- 
tion of three similar energy scales seen universally across 
members of Bi-2212 family. For the YBCO sample, we 
also can suggest the existence of the three energy scales. 
While Eo remains the same, the E\ and E2 scales are 
now shifted to lower energy relative to Bi-2212 and are 
labelled _E((not shown) and E' 2 . In addition, there is a 
possible fourth scale E' 3 ~ 0.6 eV. In the lower frame, 
we show —Y,2 P ~ qp for the four samples. This quantity 
is positive definite and shows some of the characteristic 
energy scales more unambiguously. In particular, we see 
the emergence of a new band at larger energy. In the 
case of a finite band, we expect the scattering rate to 
drop [2^, 153, HU beyond the end of the boson spectrum. 

In summary, we have demonstrated a new analysis to 
extract the quasiparticle self-energy directly from the op- 
tical conductivity data, which will facilitate future com- 
parisons with both ARPES and theory. In applying this 



method to Bi-2212 for various dopings, we find good 
agreement with the ARPES data in comparison with the 
optically-derived quasiparticle dispersion curve. While 
the optics data show features at both low and high en- 
ergies (E„ ~ 50 meV and E 2 ~ 0.9 eV) as in ARPES, 
the profile of the intermediate scale E\ ~ 0.4 eV, which 
has been identified as the background spin fluctuation 
scale, is very different in optics. Notwithstanding the 
lack of momentum resolution, optics provides a comple- 
mentary technique to ARPES in this case because it is a 
bulk probe, has enhanced energy resolution, and is able 
to examine a wider range of materials such as the YBCO 
family. Indeed, we find with the YBCO data shown here 
that we can identify up to four energy scales, which range 
from low to high energies, but that overall they appear 
to be shifted to lower energy relative to Bi-2212. The 
theoretical understanding of these features remains to be 
more fully elucidated, but confirmation of these scales 
by more than one technique gives impetus to such an 
endeavour. 
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